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Abstract 

Many models of dark matter scattering with baryons may be treated either as a 
simple contact interaction or as the exchange of a light mediator particle. We study 
an alternative, in which a continuum of light mediator states may be exchanged. This 
could arise, for instance, from coupling to a sector which is approximately conformal at 
the relevant momentum transfer scale. In the non-relativistic effective theory of dark 
matter-baryon scattering, which is useful for parametrizing direct detection signals, 
the effect of such continuum mediators is to multiply the amplitude by a function of 
the momentum transfer q, which in the simplest case is just a power law. We develop 
the basic framework and study two examples: the case where the mediator is a scalar 
operator coupling to the Higgs portal (which turns out to be highly constrained) and 
the case of an antisymmetric tensor operator that mixes with the hypercharge field 
strength and couples to dark matter tensor currents, which has an interesting viable 
parameter space. We describe the effect of such mediators on the cross sections and 
recoil energy spectra that could be observed in direct detection. 


1 Introduction 

Most of the matter in our universe, by mass, is dark matter, but beyond the fact that it 
interacts gravitationally, the nature of dark matter remains elusive. As experiments dig 
further into the parameter space of classic theories of dark matter and continue to hud null 
results, it is important that we think as broadly as possible about what dark matter might 
be and how we might detect it. In this paper, we will suggest a novel form of interaction 
between dark matter and baryons and explore the extent to which it modihes the signals 
experiments searching for dark matter might observe. 
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In recent years, a mnch wider variety of possible dark matter models and phenomenology 
has begun to be explored. Non-relativistic effective theories have systematized the explo¬ 
ration of possible operators characterizing dark matter-baryon scattering in direct detection 
experiments [l 13 , drawing on older work outside the dark matter context 14 . The basic 


operator approach can be modihed in various ways, for instance through considering dark 
matter particles inelastically scattering to or from excited states 15 17 , dark matter par¬ 
ticles with form factors 18,19 , dark matter that scatters through 2 —)■ 3 processes 20,21 


scattering of dark matter off two nucleons at once 22 -24 , or (a more radical modihca- 


tion) detection not of dark matter itself but of relativistic DM annihilation products 25 


Theories containing a large ensemble of (possibly unstable) dark matter states have been 
considered 26 -29 , as have theories in which only a small fraction of dark matter enjoys a 
richer set of interactions 30 33 . All of this theoretical exploration has helped to broaden 


our sense of what realistic theories of dark matter can be, pointing the way to new signatures 
that can be tested experimentally. 

Our goal in this paper is to explore yet another modihcation of the standard picture 
of how dark matter interacts with other particles. Specihcally, we will study dark matter 
interactions that are mediated by generic operators of arbitrary scaling dimension (consistent 
with unitarity bounds), which can be thought of as the exchange of a continuum of light 
states. Schematically, we would like to think about Lagrangians of the form 


L = yTyOmed + '0r'0Omed, 


( 1 ) 


where yTy stands in for any dark matter bilinear operator and '^T'^ for a Standard Model 
bilinear and Omed is some “mediating operator.” The precise technical meaning of the 
statement that a continuum of states is exchanged is that the momentum-space two-point 
function (C>med(<?)C>med(—<?)) has a branch cut extending down to momenta well below the 
threshold momentum exchange giving rise to a signal in the experiment, |g| <C ^exp (and 
furthermore that the spectral weight is spread out along the branch cut rather than being 
concentrated near a single narrow peak). Cases where C>med is a single light or heavy particle 
are well-studied, but the general case where it represents the exchange of a continuum of light 
states has received little attention. Such couplings to a continuum of states have appeared in 


the phenomenological literature in various guises, e.g. in the RS2 setup 34 or the literature 
on “unparticles” 35,36 , but (perhaps surprisingly) has been mostly absent from explorations 


of how dark matter can interact with the Standard Model. 

Let us mention here some related work in the literature. The possibility that direct 
detection could proceed through a continuum of states which would modify the effective 
nonrelativistic potential V{r) was briefly discussed in section 2.2 of (^, but not developed at 
length. An unusual velocity dependence of dark matter annihilation from the Sommerfeld 
effect due to exchange of a continuum of states was explored in [^; related work neglecting 
the Sommerfeld effect appeared in 38 . Finally, a related scenario involving a Randall- 


Sundrum realization of a tower of light mediators was studied in 39 43 

The outline of this paper is as follows. In Sec. we decribe our setup in detail. In 
particular, we briefly review the CFT formalism that we use in our calculations and describe 
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possible ways to model a mass gap. We show that the direct detection rates are almost 
not affected by our assumptions about the modeling of the mass gap. We also briefly ad¬ 
dress basic cosmological concerns related to our scenario. In Sec. we study continuum 
mediators coupling to the SM Higgs portal 44 48 . It is a simple illustration of the general 


idea, but we find that there is not a viable parameter space for the sort of direct detection 
signal we are interested in. In Sec. we analyze the DM direct detection and collider con¬ 
straints in the case when the DM-nucleus interaction is mediated via an antisymmetric tensor 
operator which couples to the SM via the hypercharge portal 49,50 . This case realizes in¬ 


teresting continuum-mediated phenomenology in a parameter space compatible with various 
constraints. Finally in the last section we conclude. Some technical details are relegated to 
the appendix. 


2 The scenario: mass scales and kinematics 

2.1 The basic picture 

The majority of models of dark matter-baryon scattering considered so far take one of two 
forms: 


Light particle mediator: 


Pointlike interaction: 



+ q). 


= Jxip^V - (l)JsM{,k,k + q). 


( 2 ) 

(3) 


In position space, these correspond to potentials V{r) oc 1/r and V{r) oc respec¬ 

tively. (These are point-particle idealizations and should be appropriately convolved with 
the nuclear form factor and, if it exists, dark matter form factor.) We could also consider the 
case of a massive mediator with mass m ~ g, which would correspond to a Yukawa potential 
interpolating between these two extremes. In this paper, we consider a different scenario, in 
which we exchange a continuum of light modes. One way to think of this is as the result of 
coupling to states of multiple light particles: 


Continuum mediator: 



JxiPyP-^) 



JsM{k,k + q). 


(4) 


Here the scale-invariant factor of (1/g^)" is a stand-in for more general possible behavior 
of the intermediate continuum. In the simplest case, we could consider just a loop of two 
light, non-interacting particles. More generally, the continuum could consist of multiple 
particles that are themselves interacting, as suggested by the shaded region in the figure. 
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These interactions could give the operator Omed an anomalous dimension, and in the strong- 
interaction limit could render any simple particle interpretation unreliable. 

The examples that we have written above involve a very important factorization property, 
which will apply to all of the models that we consider. Namely, the amplitude for dark 
matter-baryon scattering is a product of three factors: 


• A Standard Model current or “portal.” At the microscopic level, this may involve 
quarks or gluons. In a realistic direct detection calculation, the coupling is to nuclei, 
and so this piece of the amplitude in general involves a nuclear form factor describing 
the way that protons and neutrons are distributed within the nucleus. Such form 
factors are conveniently calculated with the code of [^. 


A dark matter current. Generally this is taken to be simpler than the Standard Model 
current, since dark matter is treated as an elementary particle. However, more gener¬ 
ally, dark matter itself could have a form factor (see e.g. 


18,51-53 


• The propagation of the mediator. In the case of a contact interaction, this factor is 
simply 1. For a massless particle, it is l/q^, corresponding to a long-range force. For 
this paper, we will take it to be of the form l/(g^)" (motivated by scale invariance) or 
more generally l/(g^ — m^)“ (which, as we will discuss below, is a crude but useful toy 
model for a mass gap). 


Importantly, the propagation of the continuum mediator that we consider will always simply 
have the effect of rescaling well-understood calculations in the literature by simple functions 
of q^. It never requires any new nuclear form factors, for instance. Furthermore, the dark 
matter and Standard Model currents may be in nontrivial representations of the Lorentz 
group, but the effect of the mediator will always resemble rescaling the result of a contact 
interaction by a Lorentz scalar function of q^. In the case where the mediator has spin and its 
propagator may involve complicated tensor structures, this is not entirely obvious. We work 
out the details for an antisymmetric tensor mediator that mixes with the electromagnetic 
held strength in Appendix showing that the amplitude is a simple product of known 
amplitudes for electric or magnetic dipole moments and a power law in —for a general 
scale-invariant mediator. 


2.2 A cosmological concern 

The simplest realization of continuum exchange is to couple dark matter and the Standard 


Model to an operator in an infrared-conformal sector 34,35 . However, this is potentially 


cosmologically dangerous, because a conformal sector behaves as dark radiation and is sub¬ 
ject to cosmological bounds due to its effect on the expansion rate of the universe 54 . These 
are usually quoted as constraints on the number of “effective neutrino species,” Ags, from 
BBN (AA"efr < 1-44 at 95% conhdence [^) and the CMB (Ags = 3.15 ± 0.23 from Planck 
combined with other data [^). If dark matter couples to dark radiation, these constraints 
can change in various ways depending on the form of the coupling 57 60 , but they remain 
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quite stringent. The safest way to avoid these bounds is if our continuum of particles does 
not extend all the way down to zero mass and develops a mass gap; if the particles become 
nonrelativistic before BBN, they are no longer dark radiation from the viewpoint of cosmo¬ 
logical bounds. For instance, our approximately conformal sector could conhne at a scale of 
at least a few MeV. This mass means that during BBN, the continuum of modes behaves as 
dark matter rather than dark radiation. We can also arrange that, after acquiring a mass, 
the continuum states simply decay before BBN, dumping their energy into the Standard 
Model plasma. For instance, if they couple to the Standard Model through dimension-six 
operators suppressed by a scale A, they can have a lifetime 


r ~ —r ~ 7 X 10 
m° 


-7 


sec 


A 


100 GeV 


^ Mo MeVy 
m J 


(5) 


easily decaying before BBN (although stable on collider timescales). More rehned estimates 
can be done for the particular models we discuss. In the particular case of the hypercharge 
portal, we will discuss such estimates further in Sec. |4.2[ It is also possible that deviations 


from thermal equilibrium just before BBN could assist in circumventing the dark radiation 
bounds even if there is no mass gap (or the mass gap is much smaller than the BBN scale) . 


2.3 The direct detection scale meets the BBN scale 


Direct detection experiments search for nuclear recoil events. If a nucleus of mass mjv recoils 
with kinetic energy Eji, the momentum transfer is g = \/2mi\rEji and the incoming dark 
matter velocity must have been at least Umin = ^ where fi is the dark matter-nucleus reduced 
mass. A typical dark matter velocity in the galactic halo is h = 220km/s, which (taking 
fi = 100 GeV for reference) can impart at most a momentum transfer of g ~ 147 MeV. A 
low-threshold nuclear recoil might be taken as, for instance, Epi = 2 keV for a silicon nucleus 
of mass ~ 26 GeV 62 , corresponding to about a 10 MeV momentum transfer. Thus, 


the range of momentum transfers that are relevant for direct detection might be roughly 
construed as 


10 MeV < g < 400 MeV, 


( 6 ) 


with the typical momentum transfer of interest in the middle of this range and the exact 
details depending on the particular experiment. LUX, for instance, studies nuclear recoils 
in xenon between about 3 keV and 25 keV 


63 , corresponding to 27 MeV < g < 78 MeV. 


From these estimates we see that, if the continuum of modes mediating scattering acquires 
a mass gap m 10 MeV in order to avoid dark radiation bounds during BBN, this mass 
will be a subdominant correction to the two-point function (Omed(?)Gmed(—?)) at the values 
of g that are most relevant for direct detection. The “mass-gap” solution to the BBN bound 
suffers from a coincidence problem: we have no explanation for why the gap should fall in 
the relatively small interval between the temperatures relevant for BBN and the momenta 
relevant for direct detection. (Particular models may offer solutions, but no general solution 
is apparent.) But if we assume that Tbbn ^ rn < gexp, we have interesting direct detection 
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phenomenology while dodging the cosmological bound. If we prefer to avoid accidental 
coincidences of mass scales, we can always pursue the more elaborate nonthermal cosmologies 
alluded to above to allow for much smaller m. 


2.4 The conformal limit 


Now that we’ve argued that we can at least approximately neglect masses over the range 
of momenta relevant for direct detection experiments, let us introduce the formalism we 
will use throughout most of the paper, which assumes the conformal held theory limit. To 
outline the basic formalism and the assumptions we rely on, we’ll discuss one of the simplest 
cases we can consider. Namely, C>med is a scalar operator of dimension d, coupling to the 
Standard Model through the Higgs portal 44,45 and to a Dirac fermion dark matter particle 
y through a scalar bilinear: 


C 


-^H^HO + 

Ad-2 


Co 

Ad-1 


Oxx- 


(7) 


(We now drop the subscript “med” for convenience.) The hrst term is marginal if d = 2; 
we could imagine generating a scalar operator with d ~ 2 in various ways, including as a 
fermion bilinear '0-0 near the edge of the conformal window in a QCD-like theory 64 . For 


most of the discussion in this paper, the only information that we need about the operator 
O is its two-point function, to which we assign a simple expression in position space. For 
future use, we also quote the result for an antisymmetric tensor operator: 


{o{x)om 


Co 


lx 


2d’ 


|x 


(a:), where 4^(x) = - 2^^. 


x^ 


( 8 ) 

( 9 ) 


Unitarity requires d > 1 for the scalar case and d > 2 for the antisymmetric tensor. In 
momentum space the two-point functions acquire extra prefactors: 


(0(-q)0(q)) ^ ld‘xe‘^-‘{O(x)O{0)) = ^^^Co(-q‘r-\ (10) 

{0^Aq)0>'’{-q)) = 

In both these expressions cq is a normalization constant and we will further take it to be 1 
to simplify the calculations. We work in a mostly-minus metric so that the branch cut arises 
at physical timelike momentum > 0. Because we will mostly use the momentum-space 
expression, one might wonder why we don’t assign it the simple coefficient cq and shift the 
unpleasant Gamma functions into the less-used position space answer. The reason is that 
we would like smooth behavior of answers near integer dimension: if d n -|- e, then 


(_^2)n 2 +e ^ ^ (l -h e log(-g^) H -) , 


( 12 ) 


6 













in which it turns out that the leading term is removed by contact terms and the log(— 
piece is physical. In fact, the log(—g^) factor is, for integer dimension operators, the source 
of the branch cut corresponding to a continuum of physical states. It appears that it would 
vanish in the e —)■ 0 limit, but this is compensated by a ^ from a pole in the prefactor 
r(2 — d), which we therefore want to keep. (Further discussion of such subtleties may be 
found in ref. [36|.) 

The next subtlety to discuss is the meaning of the scale A. We could imagine that this is 
a true cutoff: local held theory begins at the scale A, at which point operators already exist 
with unusual scaling dimensions, and we are simply given the Lagrangian. But, especially 
since the values of A that we can probe will not be far beyond the TeV scale, it seems 
more likely that the physics is established at some higher energy scale, possibly in terms 
of weakly-coupled elementary helds, and RG running toward strong coupling leads to the 
development of large anomalous dimensions and CFT-like behavior. In this case A may be 
a combination of other underlying mass scales. 

We will not dwell at length on UV completions in this paper, but let’s elaborate on this 
point. Say that the operator O has a nontrivial scaling dimension, perhaps d = 2 + e. Then 
we might parametrize its coupling to the Standard Model as -^h^'hO. This is an unusual 
expression, and especially if A is a relatively low scale (say, 10 TeV), we should expect 
that there is dynamics lurking behind it. Perhaps, for example, at some energy Aq this 
originated in a dimension-5 interaction of weakly-coupled helds, and then some 

other interactions at strong coupling drove the operator '0-0 to have a nontrivial scaling 
dimension. If these interactions became important at a scale M -C Aq, then we should 
match 'ijj'ijj —)■ and we infer that the effective scale A suppressing the operator at low 

energies is A = Aq{Aq/M)~ ^ Aq. Meanwhile, depending on the UV completion, perhaps 
there were also contact interactions directly linking the Higgs to dark matter at the scale Aq. 
The fact that A 3> Aq thus raises two concerns: 

• If M ^ Aq <C A, then we either have to make A quite large (suppressing our signals) 
or squeeze a lot of dynamics into energies near the weak scale (possibly giving rise to 
new constraints). 

• If other baryon-dark matter interactions are suppressed by Aq rather than A, they may 
give larger effects than those mediated by O. 

These concerns suggest to us that we should focus on cases where O is an operator of low 
dimension. If O mediates a long-range force, which is to say, if the effective amplitude 
for scattering mediated by O is proportional to a negative power of —g^, then the powers 
of (Ao/A) that tend to suppress the effects of (9- mediated exchange may be overcome by 
powers of the large ratio A^/(— g^) S> 1. Precisely because direct detection operates at low 
momentum transfer, the effects that we are studying may be observable. 

Given any effective continuum-mediated model, it would be interesting to follow this logic 
through more carefully in concrete UV completions. For now, however, we simply take away 
the general lesson that continuum-mediated scattering is most likely to be observable when 
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the amplitude involves negative powers of —q^, and is potentially subdominant (though this 
is UV-dependent) to contact interactions when this is not the case. 


2.5 Implications of the mass gap 


For a free held, the change from massless to massive is straightforward: we replace the 
propagator i/q^ with i/{q^ — m^). (For higher spin helds there are also modihcations to 
the tensor structure in the numerator.) But for generic operators, there is no universal 
prescription for how the mass gap appears in the two-point function. The spectral function 
will be zero below the threshold mass m, and nonzero above it, asymptoting to the gapless 
answer; but the near-threshold behavior could be quite complicated. For conhning gauge 
theories, for instance, we have no analytic tools to compute the precise spectral function. 
Nonetheless, it is useful to have some simple examples to make more quantitative statements 
about to what extent a ~ 10 MeV mass gap can alter simple scale-invariant predictions for 
direct detection rates. For direct detection, we are interested in spacelike values q^ < 0, 
because g is a momentum transfer; this helps to protect our answers from extreme sensitivity 
to the near threshold behavior, since the threshold and any associated peaks in the spectral 
function are at > 0. 

Let us explore how sensitive the answer is to different models of the mass gap. The 
simplest conceivable toy model is to replace (—g^)“ —>■ {rn? — g^)" (^. Other toy models can 
come from loops of massive particles or from conhning gauge theories or extra dimensions. 
For simplicity and concreteness, let us focus on the case of dimension 2 scalar operators, 
for which the CFT two-point function is (up to normalization) simply log(—g^). The simple 
toy model, then, would replace this with log(m^ — g^). One way to obtain a dimension 2 
scalar is as a product of two free scalars; in this case the spectral function arises from a loop 
integral. If we give the free scalar a mass mo, we introduce a branch cut for the two-particle 
operator at the threshold m = 2mo, and hnd that n(g^) oc dx log(—x(l — x)q^ + \rn?) -|- 2. 
(We shift the result by two so the asymptotics matches log(—g^) at —g^ —)■ cx).) A third 
toy model motivated by the RS2 interpretation of the conformal two-point function is to 
imagine inserting an infrared brane or “hard wall” to produce a mass gap, as in RSI [^, or 
approximately in conhning gauge theory at large ’t Hooft coupling (see e.g. [^). We can 
then compute 


n(g^) = 


iog(gV?o) - 


(13) 


If Xi ~ 2.404 is the hrst root of Jo{x), then we choose q^ = /x\ to obtain a mass threshold 

m. For this model the spectral function is a sum over poles, as in a large-gauge theory. 
In a £nite-A^ theory all of these poles acquire a width; we can approximate this effect on the 
spectral function by studying the large-answer slightly away from the real axis, effectively 
smearing out the narrow states: pa{s) = ^(n(s -|- iA) — n(s — iA)) 


69 


An alternative 

large-A model, perhaps more representative of real QCD (or other theories with small UV 
’t Hooft coupling), is the digamma function (69| - [7l] : n(g^) = V’(—g^ -|- m^), which has its 
first pole at w? and asymptotes to log(—g^) for large negative g^. 









Two-Point Functions at Spacelike Momentum Transfer Smeared Spectral Function A = 1/4 




q^lm^ 


Figure 1: Four toy models for a mass gap in the two-point function of a dimension-2 scalar 
operator. In every case, n(g^) —)■ log(—in the spacelike region ^ 0. The four models 
are similar at spacelike momentum (left) but differ greatly for timelike momentum (right). 
Solid blue: the simple replacement log(—g^-|-m^). Dashed orange: the result for a loop of free 
massive scalars. Dotted green: the Randall-Sundrum or “hard-wall” ansatz, characteristic of 
conhnement at large’t Hooft coupling. Dot-dashed red: the digamma function, a toy model 
of QCD-like conhnement. 


These four toy models are illustrated in Fig. The right-hand panel shows the smeared 
spectral function, which varies from a step function (in the simple log(—g^-|-m^) ansatz) to a 
smooth curve that turns on (the loop of massive scalars) to wildly bumpy curves characteristic 
of conhning large-iV theories with many narrow resonances. Despite these dramatically 
different spectral functions, the behavior for spacelike momentum, as shown in the left-hand 
plot, is qualitatively similar in every case. In particular, the corrections to the asymptotic 
conformal answer when —g^ S> rn? are very small. This gives us conhdence that the behavior 
of the direct detection cross section can be approximated with the simplest toy ansatz, 
g^ —)■ g^ — m?, even without a detailed model for the origin of the mass gap. Such an ansatz 
will certainly lead to the correct qualitative physics, and will be quantitatively reasonable 
unless we probe too close to |g^| ~ m?. 

Although we have only given examples for the special case d = 2, we expect that similar 
results would be obtained for other dimensions. (Hard and soft wall extra dimensional 
theories can produce results for arbitrary operator dimension by varying the bulk mass; 
generalizing the loop ansatz in a well-motivated way appears more difficult.) 


3 Higgs Portal 

First, let us briefly analyze the higgs portal for the continuum-mediated Dirac fermion DM. 
Because it involves a scalar operator, this case is formally the simplest. Although we will 
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immediately see this portal is excluded by the current data, at least within the standard 
cosmology, some of the results that we get here are relevant for a much more motivated 
hypercharge portal. The coupling of the mediator to the SM that we assume in the Higgs 
portal is 


C — ch 


Who 

A '^-2 ’ 


(14) 


with d being a scaling dimension of the scalar operator O, whose critical dimension is 1. The 
most generic coupling that one could write down to the fermionic DM is 


^ Oxx , ~ Oxx^x 

^ = coT:rT + co- 


'Ad-i 


Ad-i 


(15) 


Note that a priori the scales A in Eqs. (14) and (15) can be very different from one another. 


We absorb these differences in the couplings cq and ch which are not necessarily 0{1). 

As we claimed in Sec. these interactions can be reduced to the standard fermionic higgs 
portal with a g^-dependent coefficient. This approach is very similar to the approach that we 
are taking in the hypercharge portal scenario, where we will dehne g^-dependent moments. 
The effective higgs portal is 




with 


A 


_ chCoT{2 
^eff — 


A 


d)) 


A^-H{d) V A2 


d-2 


(16) 


(17) 


The expression for Ces is identical to that for Ces up to an obvious change Cq —)• Cq. It is 
very easy to qualitatively understand what would be the implications of this kind of higgs 
portal on direct detection. The momentum-dependent coupling Cefr(g^) would translate into 
an effective dependence of the recoil spectrum on g^, inducing an effective form factor, which 
has nothing to do with either nuclear or dark matter form factors, but rather arises because 
of the unusual properties of the mediator. 

Unfortunately, higgs portal continuum-mediated DM is not viable because of the im¬ 
mediate clash between the induced mass-gap in the mediator sector and the cosmological 


constraints. Interestingly, for any d < 4, Eq. (14) induces a mass gap in the mediator sector 
which is 


^gap — 


ChV 

Ad-2 


(4-d) 


(18) 


Of course if this mass gap is too small we can always assume extra sources, which might 
generate an additional mass gap. However, it is difficult to see how one would signihcantly 
reduce this induced mass gap without very severe hne-tuning. At Erst glance the automatic 
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presence of a mass gap is good from the point of view of satisfying BBN constraints; the 
problem is that the model is too predictive, implying too long a lifetime for a given mass 
gap. 

As we will immediately see the induced mass gap (18) is too large in those regions of 
parameter space where we can satisfy the cosmological constraints. As explained in Sec. 
the lightest particle in the mediation sector with a mass of order mgap should decay faster 
than one second. Assuming that this particle is a scalar s with a mass rris we can write down 
its matrix element as 


( 0 |C>|p) = 


d-lip-x 
s ^ 


(19) 


In this expression is an unkown 0{1) constant, that we will take to one in our further 
estimates. Therefore the term (14) induces mixing between the scalar s and the physical h, 
which the scalar s can decay through. Kinematically, only decays s —>■ e’''e“ and s —)■ 77 are 
allowed and their rates are 


r(s —)■ e^e ) 
r(s ^ 77) 


2cjrv‘^m^/ ^ 

h^e+e-{'^s) 

0^2 ^. 2 /^ 2 d —2 

_p / \ 

^2d-4^4 ^ , 


( 20 ) 

( 21 ) 


where Th^xxi^s) stands for the higgs partial decay widths if it had a mass mg, rather than 
ruh. In the relevant range of s masses, 1 MeV < uig < 100 MeV one can easily see that the 
77 channel is completely subdominant and disregard it. The decay to the electron pairs is 
the only viable decay channel of s, and because the rate is suppressed by the electron Yukawa 
squared, the constraints are not very easy to satisfy. In practice the constraint r* < 1 sec 
translates to the following bound: 


4 


A 




lOOGeVy ViOMovAa 


2d-2 


> 2 X 10 


-11 


( 22 ) 


It is easy to see that one cannot satisfy this constraint and the demand rugap < 50 MeV 
simultaneously. The latter suggests that d > 3 for any phenomenologically interesting range 
of the scale A. The former demands that d < 2.34 for A ~ 1 TeV, and even slightly smaller 
values of d for higher values of A. 

Of course, the constraint on the decay time of the lightest narrow state can be cir¬ 
cumvented if we assume non-standard cosmology, but we still view the constraints on this 
scenario as not appealing and further concentrate on the hypercharge portal, which is much 
more promising]^ Precisely because it does not involve a scalar operator, the coupling to the 
Standard Model in the hypercharge portal case does not automatically deform the CFT and 
produce a mass gap, so masses and lifetimes are no longer closely linked and there will be a 
larger viable parameter space to explore. 

^Another potential way to circumvent the constraints on the higgs portal would be speeding up the decay 
of the scalar s by introducing new couplings, which do not necessarily have anything to do with the DM, 
e.g. OF'^^. However this is not the most minimal scenario and these types of couplings come with their own 
constraints. We relegate the analysis of these non-minimal scenarios to future studies. 
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4 


Hypercharge Portals 


4.1 Portal and recoil spectrum 


We start from studying the hypercharge portal where the hypercharge held strength couples 
to an antisymmetric operator, 

^ . (23) 

The critical dimension of the antisymmetric tensor is 2, so the dimension of this operator is 
always bigger than or equal to 4. This is a unique portal in the SM, because it allows coupling 
of the BSM bosonic particles to the SM with an operator of dimension 4|^ Moreover, in the 
context of the DM, it is the only low-dimension operator which allows couplings of the DM 
to the SM which are not suppressed by the masses of heavy mediators: EW bosons, Higgs 
or BSM particles (e.g. Z'). 

The event rate and the total cross sections in the direct detection experiment strongly 
depend on the coupling of the antisymmetric tensor to the DM particles. Hereafter 
we will only consider operators with dimensions not higher than 6 in the limit d[Oij_u] —)■ 2. 
With this restriction, and assuming fermionic dark matter one can straightforwardly write 
down three different coupling of the DM to this operator: 


c = 


( 24 ) 

c = 


( 25 ) 

c = 


( 26 ) 


Because of the structure of these operators, which are identical to the magnetic dipole mo¬ 
ment (MDM), the electric dipole moment (EDM) and to the anapole by replacing —)■ 

we will loosely refer to each of these options as magnetic, electric and anapole hypercharge 
portals respectively. As in the previous sections, we do not necessarily assume that the 
couplings cb,C 2 ,C 2 and C 2 are order one. 

As we briefly outline in Sec and prove in Appendix one can reduce each of these 
couplings to the effective dipole/anapole coupling of the DM to the SM photon, reweighted 
by an appropriate power of —where the latter is the momentum exchanged between 
the DM and the nucleus in the scattering process. Therefore, each of these terms yields 
an effective (momentum-dependent) DM dipole moment, which depends on an appropriate 
power of —The explicit calculation gives the following effective magnetic dipole moment 
for the operator (24): 


/^mag 


cbC2Cos9iy r(3 — d) 

A2'^-3 4'^-2r(d + 1) 




d-2 


( 27 ) 


^Another unique portal which allows coupling of the BSM fermions to the SM with an operator of 
dimension 4 is the so-called “neutrino portal” HL. We will not consider the fermionic mediators in this 
paper, mainly because it is not easy to think about a scenario, where it would not spoil observed neutrino 
properties. 
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EDM coupling, no mass gap, /7 ?dm = 100 GeV 




EDM coupling, no mass gap, moM = 100 GeV 




Er [keV] 


Figure 2: Event rate in arbitrary units for electric-hypercharge (top panel) and magnetic- 
hypercharge (bottom panel) portals for different values of d. Green, orange and brown lines 
stand for d = 2.3, 2.7, 3.4 respectively. The curves have been artihcially scaled to intersect 
at a single point, because the goal is to convey shape information only. The isotope ^^^Xe is 
shown at left and at right. Notice that, due to different behavior of nuclear form factors, 
the shapes can be quite different for different atoms. To guide the reader’s intuition for how 
these compare to more familiar scenarios, in the upper-left hgure we also show the recoil 
spectrum arising from the standard contact operator coupling to mass (xX-NA^) as a dashed 
gray line and the same operator with a massless mediator a dotted gray line. 


For the EDM and the anapole hypercharge portals we get exactly the same expression for the 
effective moment with an obvious replacement C 2 —)■ C 2 for the EDM and C 2 —t C 2 ,A^'’*“^ —)■ 
A 2 rf -2 anapole moment. 

Because all the effective moments depend in a non-trivial way on the momentum transfer 
we expect that the differential event rates dR/dE^ will be modihed compared to the 
“standard” DM dipole/anapole scenario. One can refer to the factor as an effective 

DM form factor. However this form factor arises because of the non-trivial dynamics in 
the mediation sector, rather than non-trivial structure of the DM. Of course these effective 
DM form-factors multiply the nuclear form-factor in the direct detection picture, yielding 
distinctive event distributions in the direct detection experiments. 

It is also worth noticing that the effective DM form-factors will be further modihed 
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Xe, EDM coupling, d = 2.7 


131 




’®F, EDM coupling, cf= 2.7 




Er [keV] 


Figure 3: Event rate in arbitrary units for electric-hypercharge (top panel) and magnetic- 
hypercharge (bottom panel) portals with and without mass gap. Black (solid), blue(dashed) 
and red (dot-dashed) lines stand for rugap = 0, 15, 40 MeV respectively. We assume d = 2.7 
on all these plots. 


compared to the perfect CFT limit by the mass gap effects. If this mass gap is not too 
far away from the BBN limit, namely rugap ~ 0(1... 10) MeV, this modihcation is mostly 
important in the low-energy recoil regime. Therefore we expect stronger modihcations for 
the lighter elements, e.g. fluorine and oxygen, while the heavy elements like xenon might be 
only very weakly sensitive to the mass gap effects. 

To illustrate these points we show the expected rate distributions dR/dEn assuming scat¬ 
tering on both ^^^Xe and ^®F for both electric and magnetic dipoles on Fig. [ 2 IPI These hgures 
are meant to illustrate only the shape of the curve, and so have arbitrary normalization. 
The differences between the distributions are clearly seen. We further show the rates if the 
conformal invariance is broken by mass gaps rugap = 15 MeV and rugap = 40 MeV in Fig. 

In this case, the overall normalization remains arbitrary but the relative normalization of 
curves on the same plot is important, showing how the mass gap suppresses the overall rate 

^In xenon experiments like LUX and XenonlOO one of course is sensitive to the nuclear responses of all 
the seven stable xenon isotopes, and not only to ^^^Xe. We checked this explicitly, and the distributions that 
we show do not change significantly when all other isotopes are included. Therefore here for the illustration 
purposes we just use one isotope. 
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in a momentum-dependent manner. To model the mass gap here we use the most naive 
model, namely we replace —q^ —>■ + ?7igap)- Because the momentum transfer is spacelike 

we assume very little possible difference between different ways to model this effect (see 
Sec. 2.5 for details). As we see the effect is minor or very minor for scattering on ^^^Xe, 
however it is as expected much more pronounced for scattering on 

From here we can immediately calculate the total cross sections and the differential event 
rates in the direct detection experiments. Of course, there are differences between the electric 
dipole on one side and magnetic dipole/anapole on another side. While the cross section of 
the former is enhanced by q~^ in the d[0^i,] —)■ 2 limit, the cross sections of the anapole are 
completely regular and behave as in this limit. The MDM cross section is technically also 
enhanced by q~^ and the total cross-section diverges in the deep IR, this term is suppressed 
by the transverse velocity, translating to the non-relativistic operator ^S^-{qxv). Therefore, 
we expect to get much bigger cross sections in the EDM case than in the other cases for 
the same values of A and d. The cross sections have been calculated in multiple previous 
works on the dipole/anapole DM 


72 76 . The EDM cross section can be calculated from 


integrating the following expression 


da 


EDM 


dcosO 


/igi(g^(cos 6 ')) 

1 — cos 9 


(28) 


which is of course divergent for d = 2, but hnite for any other legitimate choice of d. 

The expression for the MDM cross section is slightly more complicated, and includes 
both a hnite dipole-dipole piece and a formally divergent dipole-charge radius interaction. 
The cross sections for these interactions are 


dcTMDM —charge 

dcos9 

C^CTMDM-MDM 

dcos9 


eVLg(?^(cos0))Z2 


dvr 

hlhLg(g^(cosg)) 


m. 


1 — cos 6 {m^ + mA)v‘^ 


(29) 

(30) 


One can get total cross sections for the DM-nucleon scattering from integrating these expres¬ 


sions over cos9. However, as we will explain later in Sec 4A, the direct detection experiments 
are sensitive to the effective cross sections, rather than “honest” total cross sections. There¬ 


fore, we will return to this question in Sec. 4.5 and get numerically the effective cross sections. 


taking into account the constraints on operator (23). 


4.2 The BBN constraint 


As we have mentioned in Section 2.2, we imagine that there is a mass gap for particles 
created by the operator and that these particles decay before BBN. First, let us suppose 
that the only particles we need to concern ourselves with are those created by Consider 
a vector particle V of mass my, which has a matrix element 


(0|C>/,^(x)|p, e) = ivrriy (e^p^, - e^p^) e 


^ip-x 


(31) 
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Here is an order-one constant whose valne can only be compnted if we nnderstand the 
detailed dynamics behind the mass gap. 

Throngh the mixing of with hypercharge, the V particle can decay into charged 
Standard Model particles. In the mass range my ~ 10 MeV that is primarily of interest 
to ns, the only kinematically accessible particles will be electrons since rrie- ~ 511 keV. 
Nentrinos also carry hypercharge, bnt we are considering a process at energies below the 
Z boson mass for which a decay to nentrinos will carry extra m\lm\ snppression in the 
amplitnde. As a result, we are interested in computing the partial width for the decay 
V —)■ e+e“. This is closely analogous to the —)■ e+e“ process in QCD, which proceeds 
via kinetic mixing of the rho and the photon, or the frequently studied case of dark photons 
(e.g. [^). The kinetic mixing of the V particle is 

6 y = 2cBey cosily . (32) 


We compute 


r(l/ -» e+e ) « ^(ycla cos^ 6wmv (^) 




(33) 


We require the lifetime t(V —>■ e^e ) < 1 s in order to avoid BBN constraints. This translates 
into a bound (dropping the phase space factors, which are negligible if my 3> 1 MeV): 



/ my \ 

UO MeVy 


> 9 X 


(34) 


As expected, then, decay of the vector can easily happen before BBN when d is not too 
large. Alternatively, we can express the bound in terms of the effective mixing as 


ey > 1.7 X 10 


-10 


10 MeV 


my 


(35) 


We note that this estimate could be overly conservative: more detailed computations of 
constraints on dark photons from BBN and the CMB have found that smaller values of ey 
may be safe 78 . On the other hand, these studies have assumed a single dark photon with 


abundance determined from thermal equilibrium; in our case, there is a whole dark sector. 
Precise constraints may depend on the full model of the sector, but we will take our simple 
estimates to be the best guide available in the absence of a detailed model-dependent study. 

The above constraint is the correct one if the only particles in the approximately scale- 
invariant sector are vectors, in which case they all decay directly through the mixing with 
hypercharge. On the other hand, familiarity with a range of examples of strongly-interacting 
theories tells us that often the lightest states are scalars or pseudoscalars. These could decay 
through other portals, like the Higgs portal—though as we noted above, the small value of 
the electron Yukawa suppresses decays through that portal. Alternatively, we could assume 
that the mixing between and hypercharge is the only coupling between this sector and 
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the Standard Model. It does not permit a direct decay of a light scalar into Standard Model 
states. However, the scale-invariant sector will in general admit conplings of snch a light 
scalar to heavier vectors, which in tnrn mix with the Standard Model. Snppose, then, that 
we have in addition to the vector V considered above a scalar S with mass ms < my- We 
expect that the strong sector will contain conplings like 


my 


(36) 


where we assnme that ms ~ my ~ mgap are all of the same order and ^svv is an order-one 


nnmber. This effective conpling allows the decay S' —)■ 77 throngh two insertions of V 
mixing. We compnte this decay width to be: 


7 


r(^ ^ 77) ~ ^svv^v cos^ GwCb ( 


my \ 

A J 


4d-8 


m. 


Anmy 


(37) 


Again, if we impose a r S 1 s constraint this leads to a bonnd: 

4rf_8 


^SVV^V^B \ ) 


ms \3 /i 5 MeV 


lO Mev) 




my J 


3.2 X lO 


-21 


(38) 


Again, we can write this as a bonnd on the mixing 


ms 


my \V2 
J Vl5MeV/ 


(39) 


In the alternative case that the lightest state is a psendoscalar P, the logic is essentially 
identical with the modihcation that the effective conpling is of the type . We take 

this to be a conservative bonnd on the constraint imposed by BBN on the hypercharge 
portal. Of conrse, if the lightest state really is a vector, we can nse the safer bonnd from 
V —>■ e~^e~. 

In both cases, we see that the most obvions danger arises when d 3 > 2. In that case the 
small ratio my/A is raised to a signihcant power and could potentially become small enough 
to cause problems. On the other hand, direct constraints tend to be strongest at small values 
of d, so we will carefully check whether the parameter space we consider below conflicts with 


these constraints. However, as we noted in Section [2721 if is also possible that a nonstandard 
cosmological history can eliminate the problem even if these bounds are violated. 


4.3 Further constraints from dark photon bounds 

Many experimental searches have placed constraints on a massive vector particle V whose 
held strength mixes with that of hypercharge. These are generally referred to as “dark pho¬ 
ton” searches and rely on a combination of precision measurements, searches for rare meson 
decays, and low-energy collider or hxed target experiments 77,79 81 . Further constraints 


on the small-mixing regime arise from the effect of dark photon production on energy loss in 
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Mixing Ey: Cb = 1/(167t^), = 1, = 10 MeV 



Figure 4: Region of parameter space allowed by the inequality (40) in blue. We have hxed 
Cb = (representative of a one-loop factor), = 1, and my = 10 MeV, and plot the 
band of allowed mixing ey as a function of d and A. 


supernovas 82 . Recent reviews of the status may be found in 83,84 . The continuum that 


we are interested in is a limiting case of a tower of dark photons with different masses and 
with a mass-dependent kinetic mixing parameter. Some constraints on dark photons carry 
over directly to our scenario. Others involve searches for narrow spectral features and no 
longer apply to a continuum of new particles. However, in the presence of a mass gap, it is 
quite plausible that near threshold there is a well-defined, narrow single-particle state. In 
the case that the spectral function consists of a sum over narrow poles from a warped extra 
dimension, it has been argued that the lightest state will give rise to the most important 


constraint 39 


One of the strongest constraints for resonances in the mass range of interest, near 10 MeV, 
is from the NA48/2 experiment at CERN, which searched for decays vr^ —)■ jV with V —)■ 
[^. They have expressed the constraint as an upper bound on ey which fluctuates, 
over the mass range 10 MeV my 20 MeV, between about 2 x 10“'^ and 8 x 10 


1-7 


Other interesting bounds come from beam dump experiments . For 


ey ~ 10 the 

Fermilab E774 experiment 88 excludes my < 10 MeV. For masses 10 MeV < my < 20 MeV, 


a combination of the SLAG E141 89 , SLAG E137 90 , KEK [9^, and Orsay 92 beam dump 
experiments exclude a wide range of mixings between about 10“® and few x 10“^. Supernova 
constraints further exclude the range 10“^° < ey ;6 10“® for dark photon masses up to 
around 100 MeV 82,93,94 and the lack of observations related to electromagnetic decays 
of the dark photon outside the supernova exclude a further range 10“^^ < ey < 10“^° for 


my < 20 MeV 95 


Reanalyzing all available data for the case of a continuum of states is both beyond the 
scope of this paper and highly model-dependent. Unlike direct detection phenomenology, 
searches for dark photons directly probe the timelike region of the two-point function for 
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and the result can be very sensitive to the detailed mechanism for generating a mass 
gap. Theories with a sequence of many narrow resonances are subject to a wider variety of 
direct constraints than theories with a broad continuum. We expect that there will likely 
be a resonance-like state near the mass threshold, based on analogy to QCD-like theories, 
but even this is not obviously guaranteed by general principles of quantum held theory. 
Furthermore, the amount of spectral weight in such a state is not guaranteed; that is to say, 
we might hnd that the parameter in (31) is signihcantly smaller than one, relaxing the 
bounds. Caveats aside, there is essentially only one window in which we might aim for the 
light state V to lie: 


2 X 10"^ < ey < 10 


-3 


(40) 


with 10 MeV < my ~ 20 MeV (to avoid the E774 constraint while still maintaining 
continuum-like direct detection phenomenology). This window squeezes in between beam 
dump constraints and the 7r° decay search. The combination of beam dump and supernova 


constraints rule out all smaller mixings above the BBN constraint (35) and are fairly robust: 


for instance, even if further dark sector states exist and a decay like V —)■ SS is allowed and 
prompt, the S particle will not promptly decay to Standard Model particles unless couplings 


through portals other than just hypercharge are added to the theory. In the window (40) of 


allowed ey, decays of the V particle are prompt on collider length scales and even the BBN 


constraint (39) for the case of scalar decays is safe. 


This constraints are visualized in Fig. We see that for reasonable values of the sup¬ 
pression scale A, in the region of allowed e (in blue) the parameter space d —>■ 2 is compatible 
with suppression scales a few orders of magnitude above the weak scale. 


4.4 Constraint from exotic Z decays 

Mixing of the Z boson with the operator allows the Z to decay into light mediators. 
These may shower or cascade in various ways, possibly producing a large multiplicity of light 
particles, like the vector V or scalar S discussed above. Depending on the lifetime of these 
particles the event could be invisible or there could be distinctive signals involving lepton 
jets 96 , photon jets 97 , displaced vertices, or combinations thereof. Without a detailed 


model, it is difficult to say precisely what the LEP constraint on such exotic decays might 
be, or even whether the tightest constraint arises from LEP rather than from a hadron 
collider. Consequently, we will focus on a very robust limit arising from the total width 
of the Z boson. The Z boson width is measured to be Pz = 2.4952 ± 0.0023 GeV 


whereas the global electroweak £t not including direct measurements of the width gives 
Tz = 2.4946 ± 0.0016 GeV [^. From these numbers we estimate the constraint on the total 
new physics contribution to the Z boson width. 


F“* < 6 MeV. 


(41) 


The bound for invisible decays is somewhat stricter, but not dramatically so 100 . 
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r(2^CFT) Constraint, ce = 1 



r(Z^CFT) Constraint: Contours of ce 
10 . 100 . 1000 . 10000 . 



A[GeV] 


Figure 5: Constraints on couplings in the theory from the Z boson exotic width. At left, we 
also £x cs = 1 and plot the largest allowed value of A as a function of the operator dimension 
d. In the region to the left of the dashed line, this value of A obeys the constraint (38) for 


BBN with a 10 MeV scalar as the lightest state; the entire plot obeys the constraint (34) for 


BBN with a 10 MeV vector lightest particle. At right, we vary A and d and plot contours of 
the largest allowed value of c^. Again, the dashed line demarcates the part of the plot that 
is safe from BBN even with a light scalar. 


We will denote the decay width induced by the mixing of the Z with by T{Z —>■ CFT), 

even though a small mass gap may mean that strictly speaking we are not decaying to a 
conformal sector. To compute this width we exploit the optical theorem, which tells us that 
F(Z —)■ CFT) = :^ImM(Z —h Z), where M(Z —hZ) is the amplitude for the Z boson to 
mix with the mediator and then mix back. The imaginary part of the amplitude comes from 
the discontinuity across the branch cut in the factor in the two-point function of 


O 




Disc = 2ism{TTd){—k‘^Y Using the identity sin(7rd) = 


27rc| SiU^ (d - 1) (d - 2) ^ 771^ ^ 

r(z ^ CFT) =- r(^)r(d + i)- I 2 A j 


2d-4 


mz- 


hnd 


(42) 


Notice that this decay width goes to zero when d —?• 2; in this limit the mixing is with a free 
particle, and the spectral function has no overlap with the Z mass. Comparing equations 
(41) and (42) gives ns a bound on the largest allowed mixing for any model. We plot this 


bound in Fig. The figure shows that at small values of d, approaching the limit of simple 
kinetic mixing, the bounds become strong and force ns to consider large values of A. At 
larger operator dimensions, the bound is weaker and even relatively low values of A are 
allowed by the data. These large operator dimensions (d > 2.9) can be in conflict with the 


BBN constraint (38) in the case that the lightest particle in the mediator sector is a scalar. 


However, as noted earlier this bound can be avoided by adding more direct couplings of the 
scalar operator to the SM (e.g. to the Higgs portal) or by unconventional thermal histories 
of the early universe. 
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4.5 Maximum total cross sections 

The constraints from dark photon searches and new contributions to the Z boson width 
lead to a maximum possible direct detection cross section for a given operator dimension 
and value of C 2 . We would like to assess how large a total cross section it is reasonable to 
obtain. However, the total cross section itself is not an ideal quantity to plot. The reason is 
that, particularly for models that produce distributions dR/dE^ that are highly peaked at 
low energy transfers, the total cross section may not be reflective of the event rate in a real 
experiment with a hnite energy threshold. (Indeed, the total cross section may not even be 
dehned, for sufficiently singular behavior at low Er.) As a result, we will define a notion of 
“effective total cross section.” 

The typical exclusion plot of an experiment like LUX gives a constraint on the dark 
matter-proton scattering cross section. Of course, this involves an assumption, and makes 
the most sense for the case in which the operator involved in the non-relativistic theory is 
simply 1, e.g. when the underlying relativistic interaction is ■ Thus, to get a rough 

notion of the total event rate allowed, we take the following definition: For a given model of 
interest, the effective dark matter—proton cross section <7^® is defined as the value of 
CTp obtained in a model with a 100 GeV Dirac dark matter particle x scattering through the 
interaction XXPP that gives the same integrated rate for scattering on ^^^Xe over the recoil 
energy range 3 keV to 25 keV as the model of interest. This definition has the disadvantage of 
making explicit reference to a particular experiment and its range of accessible recoil energies. 
However, any definition of cross section that we compare across different experiments and 
different models must build in some assumptions. Rigorously speaking, we should fit each 
individual model across all experiments. For the simple purpose of getting an order-of- 
magnitude sense of how the scattering rate allowed by our model compares to the scattering 
rate of a more standard model, on the other hand, a®® is fairly useful. 

With these interpretational caveats out of the way, we present the result for the largest 
CTp® allowed by the dark photon and Z width constraints. We compute dR/dE^ for the 

comparison model, C = j^xxPP^ using the code of Ref. 4 , then numerically integrate over 

2 / \ 2 '—* 

Eh from 3 to 25 keV. For this model, <Jp = ^ ■ We also compute dR/dE^ for our 

model of interest with the same code, putting in a contact interaction and then weighting the 
answer by appropriate powers of —q^ or of (—to model the continuum exchange. 
Finally, we scale fp so that the integrated rates match and read off the corresponding value 
of (Tp. 

The result of this exercise, for dark matter scattering on ^^^Xe, is presented in Fig. 
The plot illustrates that there is no allowed cross section when d > 2.6, because the Z width, 
beam dump, and supernova constraints forbid all possible ey ~ and smaller ey are 

excluded by BBN. For smaller values of d compatible with the range 2 x 10“^ < ey ;6 10“^, 
the largest allowed cross section is determined mostly by the dark photon constraint when 
d < 2.4 and mostly by the Z width constraint when d > 2.4. Notice that the cross sections 
allowed for EDMs are several orders of magnitude larger than for MDMs. This is because, at 
fixed coefficient cb, the non-relativistic scattering through the MDM operator is suppressed 
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MDM on ’’^'Xe: max. <7p°*[cm^], mpM = 100 GeV EDM on ’’^'Xe: max. <7p®*'[cm^], mpM = 100 GeV 
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Figure 6: The largest effective cross section a®® for dark matter scattering on ^^^Xe allowed 
by the constraints from dark photon searches and the Z boson width. The left-hand plot 
is for magnetic dipole-type couplings of DM to and the right-hand plot for electric 
dipole-type coupling. We fix C 2 = 1 and a dark matter mass of 100 GeV. For a given A 
and d, the largest allowed value of is chosen. At small d this is determined by the dark 
photon constraint ey < 10“^. At somewhat larger d 2.4 the Z width becomes the dominant 
constraint. At the point that the Z width no longer allows ey ^ 2 x 10“^, the beam dump 
and supernova constraints on dark photons force us all the way down to ey so small that the 
dark photons would decay after BBN, so no cross section is allowed above d ~ 2.6. What 
is plotted is the cross section cXp in square centimeters associated to a model with scalar 
contact interaction achieving the same integrated rate. 

by ~ 10“® relative to that through the EDM operator. 

To summarize our results, despite the existence of a variety of stringent constraints on 
the operator coupling to hypercharge, arising from both low-energy probes like beam- 
dump experiments and pion decay and high-energy probes like Z boson decays, there is a 
range of dimensions—roughly 2 < d < 2.6—for which sizable direct detection cross sections 
could occur. Of course, direct detection experiments themselves can constrain this region. 
However, due to the unusual spectral shapes that appear in continuum-mediated scattering, 
a new analysis of the direct detection data will be necessary to derive precise bounds. We 
leave such analyses for future work. 

5 Conclusions 

As the simplest models of WIMPs come under strain from a variety of experiments (direct 
detection, indirect detection, and colliders), it becomes increasingly important to broaden 
our theoretical vision of what dark matter might be. In recent years, a large number of 
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directions have been explored. One common theme is the possibility of a dark sector, with 
additional new particles beyond the dark matter alone. In this paper we have explored 
a novel type of dark sector in which the scattering of dark matter with ordinary matter 
proceeds not through a contact interaction or the exchange of a single particle but through 
a continuum of mediators. 

The basic formalism for direct detection is simple: the continuum mediator multiplies the 
amplitude by a function of which in the simplest scale-invariant scenario is just a power 
law. However, as we have seen, various complications arise. If we wish to have a mass gap 
to avoid BBN constraints, the direct detection phenomenology can be surprisingly robust, 
but the constraints imposed on the scenario from other experiments may be very sensitive 
to the nature of the mass gap. We have explored this in some detail for the case of an 
antisymmetric tensor mediator coupling to the held strength of hypercharge. A combination 
of low-energy and high-energy accelerator experiments puts signihcant restrictions on the 
allowed parameter space. Nonetheless, we have argued that these constraints still allow 
room for quite large and detectable signals at direct detection experiments like LUX. 

Many aspects of our analysis could be rehned in the future, as the interplay between 
direct detection experiments and other experiments is complex. We have not discussed 
indirect detection, where the Sommerfeld effect may be relevant (see 37 for related work). 


The annihilations are expected to proceed dominantly into the states of the approximately 
conformal sector, undergo complicated cascades and then hnally decays into the SM particles. 
Eventually the annihilation mode closely resembles the hidden valley scenario (2 —)■ many). 
This distinctive annihilation pattern together with potential signal in low-energy searches 
for the light dark photon can potentially help to differintiate this scenario from an ordinary 
form-factor DM, which otherwise looks identical to the continuum mediated DM scenario in 
the direct detection experiments. These directions would be interesting to explore further. 

Despite these complexities, the picture for direct detection is extremely simple, and 
searches for unusual dependence of signals on recoil energy are well worth pursuing. It 
is an exciting possibility that the discovery of dark matter could also be the discovery of 
a rich, interacting sector that could exhibit novel quantum held theoretic phenomena like 
scale-invariance in ways that we have not previously seen in particle physics. 
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A Factorization example: antisymmetric tensor 

Let’s consider two different scenarios that involve the dark matter tensor current operator 
coupling to some mediator. The first case is the standard electric or magnetic dipole 
moment coupling to the photon: 


Lx = + h.c. (43) 

In the second case, we couple to some more general operator of dimension d > 2 that 
mixes with the photon’s held strength; 

C 2 = CixV^'^xf^M. + h.c. + (44) 

The question is whether these two Lagrangians can lead to different tensor structures or 
form factors. We will see that they do not. 

In this computation we will use the following two propagators for the photon and the 
tensor held 0^y\ 

{Af,{q)Ay{-q)) = Yi^,y{q) = 

{0^y{q)0x„{-q)) = Pf,y,xa{q) oc 

(46) 


dfj.i' 




-[-q 


2\d-2 


dftXQua ^9lJ.X ^2 




Jvcr 9 


(45) 

- (/i z/) 


A.l Hypercharge with dipole moments 

Let’s hrst consider the case of Cx 



Because the electromagnetic held strength is dotted into the dark matter tensor current, 
we have one term where the photon propagator is 11,^^ and one where it is II^p, each mul¬ 
tiplying the appropriate momentum from the derivative from Ffj_y. A little simplihcation 
reveals that this amplitude dots the dark tensor current and visible electromagnetic current 
into the object 




(48) 
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A.2 Antisymmetric tensor mediator 

Now we consider the case of £ 2 , where the coupling is to an antisymmetric tensor held that 
kinetically mixes with hypercharge: 



In this case, we hrst have the propagator of itself, then an insertion which mixes it into 
electromagnetism, which propagates and couples to the Standard Model current. 

In this case, the structure appearing in the middle is 

. (50) 

Contracting all of the indices and simplifying, we hnd that 

= 2 (51) 

This shows that the general case of antisymmetric tensor mediator exchange is equivalent 
to the case of dark matter with dipole moments, reweighted by constant factors times an 
appropriate power of —q^ where q is the momentum exchanged between dark matter and 
the Standard Model in the scattering process. 

In particular, when the unitarity bound is saturated and d = 2, the operator is the 
held strength of an abelian gauge held and this reduces to the usual case of kinetic mixing. 
In that case, the q dependence is exactly as for ordinary dipole moment dark matter. 


A. 3 Comment 

In fact, this result should follow on general grounds. The amplitude necessarily has the 
form where is antisymmetric in fi and z/. The only possible tensor 

structure is p(g). More generally, exchange of a mediator held will always produce at most 
a small hnite set of tensor structures coupling a dark current to a Standard Model current. 
If an analysis of the full set of tensor structures has been carried out for the exchange of 
ordinary weakly-coupled particles, then the more general exchange of operators of arbitrary 
dimension will not lead to new tensor structures or form factors, but only to new momentum 
dependence in the amplitude. 
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